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ABSTRACT 

The local power law relationship between the surface densities of neutral hydrogen gas and star formation rate 
(SFR) can be used to explore the SFR properties of damped Lyman a (DLA) systems at higher redshift. We 
find that while the SFR densities for DLA systems are consistent with luminous star forming galaxies at redshifts 
below z w 0.6, at higher redshifts their SFR density is too low for them to provide a significant contribution to 
the cosmic star formation history (SFH). This suggests that the majority of DLAs may be a distinct population 
from the Lyman break galaxies (LBGs) or submillimeter star-forming galaxies that together dominate the SFR 
density at high redshift. It is also possible that the DLAs do not trace the bulk of the neutral gas at high redshift. 
The metallicity properties of DLAs are consistent with this interpretation. The DLAs show a metal mass density 
lower by two orders of magnitude at all redshifts than that inferred from the SFH of the universe. These results are 
consistent with DLAs being dominated by low mass systems having low SFRs or a late onset of star formation, 
similar to the star formation histories of dwarf galaxies in the local universe. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: starburst — quasars: absorption lines 
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1. INTRODUCTION 

The nature of damped Lya (DLA) galaxies remains an un- 
resolved question, although recent work has advanced our un- 
derstanding considerably. York et al. (1986) proposed that at 
least some DLAs were dwarf and low surface-brightness galax- 
ies. This has been confirmed in a number of cases for DLAs at 
z < 1 (e.g., Schulte-Ladbecket al. 2004; Rao et al. 2003; Turn- 
shek et al. 2001; Rao & Turnshek 1998). There are alternative 
suggestions that DLAs are similar to the population of Lyman 
break galaxies (LBGs) at high redshift based on their star for- 
mation (SF) properties (Wolfe et al. 2003a,b; Weatherley et al. 
2005). Here we explore the nature of the DLA population fur- 
ther through their star formation history (SFH). Using the local 
relation between gas and star formation rate (SFR) surface den- 
sities from Kennicutt (1998) we compare the space densities of 
DLA SFRs with the cosmic history of SFR density for luminous 
galaxies. 

The evolution of the volume density of SFR in galaxies has 
been recently summarised by Hopkins (2004, his Figure 1, for 
example), and we reproduce those data points here in Figure 1 . 
The hatched envelope in this figure encompasses the majority 
of these data, including estimates corresponding to the "SFR- 
dependent" obscuration correction, which increases slightly the 
height of the upper envelope above z w 1. This envelope can 
be used as a realistic bound from observational measurements 
to constrain the cosmic star formation history (SFH) and its re- 
lated integral quantities, and its vertices are given in Table 1. 
The parameterisation of the SFH (corrected for dust obscura- 
tion) from Cole et al. (2001) has come to be relied upon by the 
community due to its useful analytic form, and it is illustrated 
in Figure 1, shown as the dashed curve. This curve, based on 
the much smaller data compilation available at that time, is con- 
sistent with the current compilation at z < 2 although it signifi- 
cantly underestimates the SFH for higher redshifts. Also shown 
as the solid line in Figure 1 is an updated fit using the Cole et 
al. (2001) analytic form, = (a + bz)h/(l + (z/c) d ), with 
h = 0.7, and parameters a = 0.02, b = 0.14, c = 3.6 and 
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d = 3.4. 

A Salpeter (1955) initial mass function (IMF) is assumed 
throughout, with lower and upper mass limits of 0.1 and 
100 M©. We adopt a cosmology with h — 0.7, ft^j = 
0.3, Q A = 0.7, where H = 100 /ikms" 1 Mpc" 1 . 

2. SFR DENSITY OF DLAS 

Kennicutt (1998), using Ha, Hi, CO and FIR measurements 
for a sample of local spiral and starburst galaxies, quantifies 
a global Schmidt law relating the disk-averaged surface densi- 
ties of gas and SFR. This relationship, Equation 4 of Kennicutt 
(1998), can be expressed as: 
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Since the gas surface density in DLA systems is just the column 
density, this relationship, assuming it is valid at the redshifts 
of the DLA absorbers (as suggested by Lanzetta et al. 2002), 
can be applied to estimate the SFR in the DLA systems. First 
consider the Hi mass density pm = (Ho/c)mn J Nf(N)dN, 
where is the Hi column density and f(N) is the Hi column 
density distribution. This quantity is most commonly used in 
the calculation of Qdla'- 
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Here tor is the mass of a hydrogen atom, p CI ;t is the critical 
density, and p — 1 .3 is the factor commonly used to incorporate 
helium as well as hydrogen in the estimate of the DLA neutral 
gas density (e.g., Rao & Turnshek 2000). To calculate the SFR 
density, then, it is necessary to integrate the SFR surface density 
(derived from the column density using Equation 1) multiplied 
by the number of systems f(N) of the corresponding column 
density, over N: 
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This gives p* in units of M Q yr -1 Mpc~ 3 for mn in units of 
M Q , N in units of Mpc~ 2 and f(N) in units of Mpc 2 . For our 
assumed cosmology, this can be rewritten as 

,«„ 

p* = 4.29 x 1(T 30 / N 1A f(N)dN (4) 

•/iVmin 

with N and f(N) in traditional units of cm -2 and cm 2 , re- 
spectively. The column density distribution function is often 
parameterised as f(N) — BN@ , giving 

P* = 4.29 x 10 - 30 ^ (C - • (5) 

The quantity f(N) requires detection of many absorber sys- 
tems to be reliably measured, and is thus difficult to measure as 
a function of redshift. Conveniently, though, 
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where noLA is the redshift distribution of the DLA systems 
dn/dz, and the "absorption distance" dX/dz = (c/Hq)(1 + 
zf/E{z), with E(z) = (SI M {1 + zf + n k (l + z) 2 + fi A ) - 5 . 
This allows the mean column density (Nm) to be explicitly 
used, rather than the integral of f(N), when estimating mass 
densities of DLA gas: 

™DLA ,, r , 

Phi = j^rnr mu{N m ). (7) 



dX/dz 
We could analogously write 
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Alternatively, from Equation 8 we could invoke Equation 1 to 
infer 
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(specifying that mu(Nui) is in units of M Q Mpc -2 ). Here we 
are making the coarse assumption that the mean column den- 
sity can be used to infer the mean SFR surface density for the 
whole population of DLAs. We will show that this assumption 
is not strictly true, but it does provide a useful approximation, 
and allows estimates of p* to be made in the absence of robustly 
measured f(N) distributions. 

In order to evaluate the reliability of Equation 10, we need to 
show that the ratio 
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From Equations 6 and 9 we have 
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Now, with values of B, (3 for f(N) distributions at given red- 
shifts, along with tidla? we can evaluate R corresponding to in- 
tegrals over appropriate ranges of N. At intermediate redshift, 
0-1 < z < 1.6, f3 = -1.40, B = 10 6 ' 40 , and n DL A ~ 0.1 (Rao 
et al. 2005b), integrating over the range 20.3 < log (AT) < 21.8 
gives R = 1.4. At high redshift, 2 < z < 4, f3 = -1.79, 
B = 10 14 ' 8 , and n D LA ~ 0.25, and integrating over the same 
range gives R = 1.3. At z = 0, f(N) has been determined by 
Ryan- Weber et al. (2003) and parameterised as a broken power- 
law, with (3 = -1.4 for 19.6 < log(iV) < 20.9 and (5 = -2.1 
for 20.9 < \og{N) < 21.6. The corresponding normalisations 
are B = 10 6 08 and B = 10 2146 . Locally, n DL A = 0.046 
(Ryan- Weber et al. 2005). Then integrating over the full range 
19.6 < log(iV) < 21.6 gives R= 1.6. 

The advantage of Equation 10 is that it allows published val- 
ues of SIdla and wdla to be used to derive p* . This also allows 
a higher redshift resolution to be obtained for p*, as ^dla and 
nDLA have been estimated for more redshift bins than the num- 
ber for which f(N) has been accurately parameterised. For an 
arbitrary cosmology: 
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where cl and c2 represent the measurements using the initial 
and desired cosmological parameters, respectively. The term in 
the brackets is just toh(-^hi) derived from the published pa- 
rameters. Where no cosmological conversion is necessary, this 
simplifies to 
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We find that Equation 15 actually gives different estimates 
for p* than the equivalent calculation from explicitly using 
Equation 6 prior to applying Equation 10. This is in the sense 
that Equation 15 produces larger values of p*, and we explain 
this by noting that the observed f(N) measurements (used to 
estimate ^dla) deviate somewhat from a simple power law. 
The p* estimates from Equation 15 are coincidentally very sim- 
ilar to the estimates inferred from explicitly evaluating Equa- 
tion 3. What this means is that if the true f(N) distribution 
really is a power law, and this is well-modelled by the assumed 
parameters detailed above, then the results inferred from Equa- 
tion 15 are appropriate estimates of the true p* values. If the 
f(N) distributions instead do have significant deviations from 
a power law, then the values for p* from Equation 15 will need 
to be scaled by the appropriate R value above to reliably ap- 
proximate p*. 

Bearing this in mind, we obtain estimates for p* at moder- 
ate to high redshift using Equation 15, with values of SIdla 
and udla taken from Rao et al. (2005b) and Prochaska & 
Herbert-Fort (2004). At z = 0, where f(N) is more reliably 
determined, we use Equation 3 with f(N) as parameterised by 
Ryan- Weber et al. (2003). These are shown in Figure 2 and 
compared with the envelope defined in §1, corresponding to p* 
estimates from luminous galaxies. This result is explored in 
detail in 55 below. 
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3. STELLAR MASS DENSITY 

The time integral of p\ establishes the stellar mass density, 
p*(z), as a function of redshift. The envelope showing this in 
Figure 3 is derived from the hatched region in Figure 1, and a 
mass lock-up fraction of 0.72 has been used (Cole et al. 2001), 
appropriate for the Salpeter (1955) IMF. In this diagram, we 
assume the redshift of formation is Zf = 10, although the re- 
sults are not strongly sensitive to this choice. The Cole et al. 
(2001) parameterisation, along with our new fit, are similarly 
integrated. The differences between the Cole et al. (2001) curve 
and the high-z measurements in Figure 1 appear here as distinct 
predictions of the stellar mass density for z > 1. 

Several independent measurements of p*(z) are shown for 
comparison, after conversion to our adopted cosmology where 
necessary. At low redshift (z < 1), the measurements are en- 
couragingly consistent with the integral of the SFH (Brinch- 
mann & Ellis 2000; Cole et al. 2001; Rudnick et al. 2003), 
albeit toward the lower limit of the envelope. The data from 
Brinchmann & Ellis (2000) have been increased by a factor 
1.25, as in earlier comparisons (Dickinson et al. 2003; Rudnick 
et al. 2003), to account for their claimed « 80% completeness. 
Above z « 1.5, there is a clear discrepancy between the SFH 
integral and the measurements from Rudnick et al. (2003) and 
Dickinson et al. (2003). It is important to emphasise that the 
Salpeter (1955) IMF has been used in all measurements shown 
here, since the resulting values are strongly sensitive to the 
assumed IMF. 

The data of Dickinson et al. (2003) are shown with uncertain- 
ties representing the larger of their 1 a or systematic uncertain- 
ties for each point. Their mass estimates are based on SED tem- 
plate fitting to a Hubble Deep Field (HDF) photometric sample, 
with upper limits derived by invoking a two-component model 
that includes a maximally old, but unobscured stellar popula- 
tion. They comment that arbitrarily large masses could in prin- 
ciple be derived by freely adding extinction to this component. 
Although this aspect is not further pursued by those authors, 
this possibility may not be unreasonable. Indeed Dunne et al. 
(2003) have shown, from a dust-mass function estimate based 
on an analysis of submillimeter sources, that there is significant 
mass density in dust at z « 2.5, at least comparable to the mass 
density in stellar objects at that epoch. This suggests that heav- 
ily obscured old stellar populations might well be significant, 
and it is not unreasonable that the estimates from the HDF may 
be low by at least a factor of two. This would bring the upper 
limits of their uncertainties to within the lower region of the 
SFH integral envelope. The FIRES measurements of (Rudnick 
et al. 2003) are scaled to account for incomplete sampling, as- 
suming that the ratio between their SDSS derived estimate at 
z = 0.1 and the measurement of Cole et al. (2001) applies to 
the higher redshift data. This type of correction is highly sen- 
sitive to assumptions about the evolution of the faint end of the 
luminosity function, and it is again not unreasonable that the 
total mass has been underestimated. 

This discrepancy between the observed values for p*(z) and 
the integral of the SFH has been explored in some detail by 
Nagamine et al. (2004a). They find that simulations and theo- 
retical models produce a stellar mass density at z « 3 higher 
than observations, but consistent with the SFH. They interpret 
this to suggest that observations might be missing almost half 
the stellar mass at high redshift, citing incomplete galaxy pop- 
ulation sampling, and/or cosmic variance affecting the surveys 
that examine only small fields of view. Encouragingly, the evo- 



lution of (z) found here from the integral of the observed 
SFH is consistent with that from the simulations described by 
Nagamine et al. (2004a). It is important to verify that the enve- 
lope derived from the integral of the SFH is consistent with in- 
dependently measured data for p*(z), for two reasons. First, to 
ensure that the SFH envelope itself is a robust constraint on the 
true SFH, and second because we are about to employ it to ex- 
plore the metal mass density as a function of redshift. The latter 
reason is worth emphasising particularly, as there are relatively 
few independent measurements of this quantity at high redshift 
to serve as consistency checks on this aspect of the analysis. 

4. METAL MASS DENSITY 

For some time it has been recognized that the low mean cos- 
mic metallicities observed in DLAs are inconsistent with ob- 
served stellar metallicities. For example, initial recognition 
of this so-called "missing metals" problem led Lanzetta et al. 
(1995) to discuss it in terms of the "cosmic G-dwarf problem" 
in the context of a "closed-box" scenario for galaxy formation. 
In particular, it was clear that if DLAs traced only luminous 
disk galaxies, their mean cosmic metallicities at moderate to 
high redshift should be nearly an order of magnitude higher 
than the observed values of w 0.1 solar. Higher observed 
metallicities would be required to match the typical metallic- 
ities of solar-type stars in the Milky Way disk since they were 
formed 5 Gyrs ago. As a solution to this problem, Lanzetta et 
al. (1995) proposed that DLAs trace not only the evolution of 
galactic disks but also the evolution of galactic spheroids. In 
part this conclusion was motivated by the supposed constraint 
that the cosmological DLA gas mass density at z = 3.5 was 
converted to the mass density of stars at z = 0; but this con- 
straint is no longer required by the data. More recent sugges- 
tions that also lead to the expectation that DLA metallicities 
should be lower include the possibilities of significant neutral 
gas in dwarf galaxies and low surface brightness galaxies (e.g. 
Rao et al. 2003) as well as significant neutral gas at large galac- 
tocentric distances coupled with reasonable metallicity gradi- 
ents (Chen et al. 2005). At the time of the Lanzetta et al. (1995) 
study, however, there were no measurements of the cosmic SFH 
to further constrain the problem. Now, with better measure- 
ments of the stellar and gaseous components as a function of 
redshift (§2), the nature of the missing metals problem can be 
defined more clearly. 

The problem has also recently been explored by Dunne et 
al. (2003) through the partitioning of metallicity in the various 
components of the Universe at z = and z = 2.5. They con- 
clude that the DLAs at high redshift are not the same population 
as the dusty, highly star-forming submillimeter galaxies, which 
contain the majority of the metals at this epoch (and possibly 
earlier). 

In the context of the current study, the metal mass density, 
pz, as a function of redshift can be established from the SFH, 
since is related to pz (e.g., Madau et al. 1996). More recent 
stellar population synthesis results (Bruzual & Chariot 2003) 
indicate this relationship is p\ = 63.7pz (see, e.g., Conti et al. 
2003), and the pz so derived is shown in Figure 4. The local 
value of p z can be compared with the compilation of Calura & 
Matteucci (2004), who favor a value of 1.31 x 10 7 M Mpc -3 , 
toward the low end of the range. Values at z = and z — 2.5 
from Dunne et al. (2003) are also shown, indicating that the 
evolution in pz from the SFH is reasonably consistent with that 
estimated from the dusty submillimeter galaxy population. The 
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mass density in metals inferred from the DLAs, also shown 
in Figure 4, give a significantly different result in compari- 
son to the estimate from the SFH. To determine the mass den- 
sity in metals we use the metallicity measurements from Rao 
et al. (2005a) relative to a solar metallicity mass fraction of 
Z Q = 0.02 and log(ft D LA) = -3 for z > 0.6 from Rao 
et al. (2005b), which corresponds to log(rini) = -3.11 or 
Phi = 1-05 x 10 8 M Q Mpc- 3 . Thus, since Rao et al. (2005a) 
find that the mean cosmic metallicity of DLAs evolves from 
[M/H]= -1.6 at z = 4 to [M/H]= -0.9 at z = 1, we infer 
that the mean cosmic metallicity measured directly from DLAs 
at z > 0.6 is about two orders of magnitude less than what 
is inferred from the SFH. Figure 4 clearly illustrates the cur- 
rent status of the so-called "missing metals" problem in a way 
that does not easily lend itself to solutions involving spheroids, 
dwarfs, low surface brightness galaxies, or gas at large galacto- 
centric distance. 

Thus, based on the comparison presented in Figure 4, we 
suggest that the missing metals problem may in fact be the re- 
sult of missing a substantial fraction of the metal-enriched gas 
in DLA surveys. This gas might be in either neutral or molec- 
ular form, but even if it were in molecular form a DLA-size 
column of gas would be intercepted. In particular, the method 
of identifying DLA galaxies is based on gas-cross-section se- 
lection, and one possibility is that the neutral or molecular re- 
gions containing most of the metals have very small gas cross 
sections, leading to a situation where we are missing significant 
amounts of both gas and metals. In §5 we discuss this possibil- 
ity further. 

5. DISCUSSION 

Is it valid to apply the local relationship between SFR and gas 
surface densities to DLA systems at high redshift, as suggested 
by Lanzetta et al. (2002)? It seems reasonable that the physical 
basis underlying the Schmidt law (Kennicutt 1998, and refer- 
ences therein) would not change significantly with time, and 
the broad range of gas mass densities over which the relation- 
ship has been characterized locally support this idea. But DLA 
systems may be comprised of quite different gas components 
compared to local star forming galaxies. There are suggestions 
that very little molecular gas resides in detected DLAs (e.g., 
Curran et al. 2004a,b). Such a lack of molecular gas would 
suggest that if the Kennicutt relation does not hold for DLAs, 
it is likely that the overall SFRs would be even lower than that 
inferred here. 

There is, moreover, a discrepancy between the SFH deter- 
mined here and that inferred by Wolfe et al. (2003b) using 
the Cll* technique. Even if the R = 1.3 scaling (§2) is ap- 
plied to our p\ estimates at z > 2, this is an increase of only 
0.11 dex, not enough to significantly reduce the discrepancy. In 
the currently adopted cosmology their "consensus" model gives 
p* = 0.21±g;^ at z = 2.15, and p* = 0.19±g;n at z = 3.70 
(Figure 2). These high values are marginally consistent with 
the current result, but are more consistent with the SFH of lu- 
minous galaxies. The results are interpreted by those authors 
as evidence that the DLA population is the same as the LBGs. 
What could cause such a significant difference in derived /?*? 
There are a number of possible reasons, some of which may be 
related to the assumptions underlying the C II* SFR estimates. 
These include the requirement of a steady state condition (cool- 
ing from C II* A = 1335.7 absorption used to infer the heating 
rate), as well as assumptions regarding the abundance ratios for 



[Fe/Si] and [Si/H], and grain composition and depletion ratio 
details for different dust models (Wolfe et al. 2003a). All of 
these independently may be reasonable approximations, but to- 
gether it is possible that the combination of individual uncer- 
tainties introduces an overall uncertainty larger than that ac- 
counted for in the formal error analysis. There are also possible 
issues arising from observational biases and effects, if the ab- 
sorption lines are marginally saturated, for example, and from 
the inevitable bias toward measuring the largest contribution to 
the DLA gas cross section. Conversely, the p* estimates from 
the Hi column densities rely on the single, reasonable assump- 
tion that the Kennicutt relation is valid for the DLA systems. 
We note that the recent estimate of the external ultraviolet (UV) 
background radiation (Miralda-Escude 2005) is higher than as- 
sumed by Wolfe et al. (2003a), and would lead to lower SFR 
densities inferred from the C II* method. In a more recent anal- 
ysis (Wolfe 2005), using a UV background estimate consistent 
with that of Miralda-Escude (2005), the C II* estimates of p* 
are indeed found to be lower than those of Wolfe et al. (2003b), 
and now appear highly consistent with the results derived here. 

We turn now to related properties of the DLA population, in 
order to try to explain their apparently low average SFRs at high 
redshift. The mass densities in stars, p*, and DLA neutral gas, 
Pdla, are compared in Figure 5, along with the total baryonic 
density for reference. It can be seen that the gas currently de- 
tected in DLA systems evolves very little over the entire range 
< z < 5, even including the factor of <~ 2 evolution between 
z = and z = 0.6. Also shown in this Figure is a simple esti- 
mate of how the neutral gas mass density, p gas , in galaxies may 
evolve (see also the closed box model of Pei & Fall 1995). This 
is obtained by assuming that the total p* (z) + p gas (z) in galax- 
ies at all epochs is the same as the local value, p* (0) + p gas (0) 
(see also Salucci & Persic 1999, their equation 5). We thus have 
Psas(z) = p*(0) + Pgas(O) - p*{z). Interestingly, the density 
of neutral gas in DLA systems under this assumption appears 
to be a small fraction of the total p gas in galaxies at redshifts 
above z m 1. This difference is exacerbated if the stellar den- 
sity evolution from Cole et al. (2001) is assumed. This picture is 
quite different from the widely accepted assumption that DLA 
systems host the majority of the neutral gas at all redshifts. Al- 
though some recent work suggested that sub-DLA systems may 
contribute almost 50% to pdla at z > 3.5 (Peroux et al. 2003), 
this now appears not to be the case (Prochaska & Herbert-Fort 
2004). 

How else can we explain this apparent discrepancy? It is un- 
likely that the difference can be attributed to DLAs missed due 
to dust obscuration in quasars, since the limit on the contribu- 
tion from missed DLAs is at most around a factor of two (Elli- 
son et al. 2001, 2004). One possibility suggested in §4 is that 
there is significant mass in very high column density regions, 
corresponding to the high SFR regions in LBG or submillimeter 
galaxy systems, but having such small physical cross-sectional 
areas that they are too rare to have been detected in DLA sur- 
veys, which have so far led to the discovery of a few hundred 
DLAs (and about 100 have metallicity measurements). This 
is consistent with the SFR intensity distribution function illus- 
trated by Lanzetta et al. (2002), which suggests that DLA sys- 
tems only contribute to the low SFR surface density end. In 
turn this implies that the high SFR surface density objects cor- 
respond to very high column density neutral plus molecular gas 
(« 10 22 to « 10 25 atoms cm -2 ), which is up to three orders of 
magnitude higher than the highest measured DLA column den- 
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sities. These column densities are consistent with the high gas 
and SFR surface densities seen in local starburst nuclei (Kenni- 
cutt 1998). 

We can make a rough estimate of the requirements which 
would lead to us missing significant metal-rich neutral and 
molecular gas in current DLA surveys. The data suggest that 
we need to find w 10 times more gas, and that it needs to have 
a metallicity > 0.1 times solar near z w 4 and greater than 
solar metallicity near z « 1. At the same time, since sev- 
eral hundred DLAs have been identified (but only about 100 
have measured metallicities), this rare metal-rich population of 
DLA absorbers should have an incidence which is < 1% of the 
known DLA incidence. In current DLA surveys the mean neu- 
tral hydrogen column density of a DLA is « 10 21 atoms cm -2 
(Rao et al. 2005b). Therefore, the required population of metal- 
rich DLAs should have (iVni) > 10 24 ro.oi atoms cm -2 where 
ro. oi = 1 is the current conservative upper limit on the ra- 
tio of the incidence of metal-rich DLAs to DLAs relative to 
1%. We note that a single metal-rich absorber with an effec- 
tive radius of « 100 pc (e.g., typical of giant neutral hydrogen 
clouds in galaxies and a factor of a few larger than typical giant 
molecular clouds which are the sites for SF) has a cross sec- 
tion which is w 10 4 times smaller than known DLAs, which 
typically have effective radii of « lOkpc. Thus, if r .oi = 1, 
only with a survey that identified w 10 3 DLAs would we detect 
« 10 absorbers from this putative metal-rich population with 
(A^hi+h 2 ) ~ 10 24 atomscm~ 2 . The rare z ss 0.6 absorber in 
Q02 18+357 (Muller 2005) may be such an example. 

Another alternative, of course, is that the closed box assump- 
tion of a fixed p*(z) + p gas (z) in galaxies at all redshifts may 
be too simplistic. The complexities of galaxy evolution include 
continuous gas infall to galaxy potential wells, and subsequent 
cooling to convert high redshift ionized gas in the IGM to lower 
redshift neutral gas in galaxies. This mechanism may contribute 
to the relative constancy of /?dla, with the neutral gas being re- 
plenished as star formation progresses in these systems. Such a 
scenario, however, might very well still lead to inconsistencies 
with observed stellar metallicities. 

In any case, the low value of pdla compared to the stellar 
mass density in luminous galaxies is one of the primary rea- 
sons for the disparity in pz shown in Figure 4. The combination 
of low average metallicity plus low gas mass density naturally 
leads to a space density of metals significantly lower than in the 
luminous galaxy population. 

Finally, consider the properties of the observed population 
of DLAs in isolation. Within the hierarchical scenario, results 
on DLA evolution might be explained by assuming the popula- 
tion is dominated by dwarf galaxies, especially at high redshift. 
Dwarf systems have low SFRs on average (e.g., Hopkins et al. 
2002), or a late onset of SF (Hopkins et al. 2001; Mateo 1998), 
as well as column densities within the DLA regime (e.g., the 
Small Magellanic Cloud, Tumlinson et al. 2002). They also 
have low metallicities which are maintained at a low level de- 
spite SF events through galactic winds or other outflow (e.g., 
Mac Low & Ferrara 1999). The DLAs appear to contribute 
significantly to the SFH only for z < 0.6, which would corre- 
spond to the epoch at which the late onset SF would occur. At 
least some moderately low redshift DLAs (z w 0.5) are likely 
to be near L* disk galaxies (Rao et al. 2003), and this is also 
consistent with the DLA systems appearing to contribute more 
significantly to the SFH around this redshift. 

It is also possible that DLAs, whether dwarfs or not, could 



exist in the outskirts of the same halos as the more centrally 
located LBGs. This would naturally account for the column 
density of DLAs being lower than the LBGs, and would also 
result in a higher interception probability for the DLAs. The 
low metallicity of the DLA systems could also be naturally ex- 
plained as a result of metallicity gradients in the halo (Chen et 
al. 2005). This would further account for the DLAs not dom- 
inating the neutral gas content, because most of the gas would 
be contained in the central part of the halo with the LBGs. This 
scenario then makes a strong prediction that LBGs should be 
seen in the vicinity of all or most DLAs. The relative scarcity 
of such alignments (only one has been reported in the literature 
to date, Muller 2005, despite significant efforts to image DLA 
galaxies) suggests that this scenario is unlikely to be represen- 
tative of the majority of DLA systems. Perhaps a more rea- 
sonable scenario is that outlined by Mo et al. (1998, 1999), in 
which DLA systems are the progenitors of disk and low-surface 
brightness galaxies, evolving in predominantly low galaxy den- 
sity environments (Mo et al. 1998), while LBGs are the pro- 
genitors of massive early type galaxies, evolving in high galaxy 
density environments (Mo et al. 1999). 

The nature of DLA systems has also been explored through 
simulations. Results from Okoshi & Nagashima (2005) and 
Nagamine et al. (2004b), for example, indicate that the average 
masses for these systems are likely to be low, perhaps around 
10 9 Mq, which is consistent with a galaxy population having 
low integrated levels of SFR. The SFR distribution for DLAs 
was explored by Okoshi & Nagashima (2005) who find a broad 
range, spanning 10~ 6 to 100M Q yr _1 , but with a mean value 
of 0.01 M Q yr _1 for DLAs lying within < z < 1. They con- 
clude that DLAs are dwarf systems with typically low SFRs. 
This is consistent with the SFH seen in the current study, and 
also suggests a plausible reason for the low pz seen in the DLA 
population. Given the outflow mechanism for maintaining the 
low metallicity in dwarfs, the stellar populations in the DLA 
systems (as well as their ISM) would be expected to show the 
low metallicities characteristic of nearby dwarf systems. This 
is indeed seen in a local (z s=s 0.01) dwarf galaxy identified 
as a DLA (Schulte-Ladbeck et al. 2004, 2005). These results 
are not inconsistent with the high SFRs estimated from [OIII] 
emission in three high-z DLAs recently detected by Weather- 
ley et al. (2005). These systems are likely to be at the bright 
end of the distribution of DLA properties, consistent with the 
high SFR tail seen in the simulations of Okoshi & Nagashima 
(2005). This is further suggested by a search for Ha emission 
from DLAs at z > 2, in which Bunker et al. (1999) find upper 
limits on the SFRs of a few tens of solar masses per year. 

6. SUMMARY 

Using the local relationship between gas and SFR surface 
densities from Kennicutt (1998) we derive the SFH of the DLA 
population. At low redshift (z < 0.6) the DLAs appear to con- 
tribute significantly to the SFH, but they have a much lower 
contribution at high redshift. This indicates that at high red- 
shift the majority of the DLA systems are unlikely to be the 
same population as the LBGs, or at least they sample a very 
different luminosity regime. An exploration of the evolution 
of stellar and metal mass densities, as well as a comparison of 
Pdla with the inferred total gas mass density, suggests that the 
DLA population may be dominated by dwarf like systems with 
low average SFRs or a late onset of SF. It also suggests that the 
DLAs identified so far may not account for the majority of the 
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neutral gas at high redshift. 
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Table 1 



Vertices of the envelope in Figure 1. 



Vertex 


log(l + z) 


z 


log(p*) 


A 


0.00 


0.00 


-2.02 


B 


0.48 


2.02 


-0.90 


C 


0.68 


3.79 


-0.90 


D 


0.83 


5.76 


-1.05 


E 


0.83 


5.76 


-0.35 


F 


0.36 


1.29 


-0.35 


G 


0.00 


0.00 


-1.52 



Note. — Vertices are labelled 
ABCDEFG moving anticlockwise about 
the envelope, starting from the lower z = 
vertex. 
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log(l + z) 

FIG. 1 . — The compilation of SFR density data from Hopkins (2004). The dashed line is the parameterised form from Cole et al. (2001), which is consistent with 
the data only for z < 2. At higher redshifts this curve underestimates the SFR density compared with the measurements. The solid line is an updated fit to the data 
using the analytic form from Cole et al. (2001) as detailed in the text. The hatched region, established by visual inspection, encompasses the majority of the most 
reliable SFR density estimates at all redshifts. 
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FIG. 2. — The SFR density of DLA systems (filled circles) as a function of redshift, derived as described in the text. The crosses are the estimates from Wolfe et 
al. (2003b) using the C II* method. Solid line and hatched region as in Figure 1. 
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z 

FIG. 3. — The hatched region shows the mass density in stars derived from integrating the SFH corresponding to the hatched region in Figure 1. The dashed and 
solid lines come from integrating the SFH of Cole et al. (2001), and the new fit using that form, respectively. The open circle is the local stellar density from Cole et 
al. (2001); the filled circle and filled squares represent the SDSS and FIRES data, respectively, from Rudnick et al. (2003), scaled such that the SDSS measurement 
is consistent with that from Cole et al. (2001); the open stars are from Brinchmann & Ellis (2000); and the open squares are from Dickinson et al. (2003). 
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z 

FIG. 4. — The hatched region shows the mass density in metals derived from the SFH corresponding to the hatched region in Figure 1. The filled circles are the 
DLA measurements, from Rao et al. (2005a); the triangles at z = are the data from Table 7 of Calura & Matteucci (2004); the open circles are from Table 1 of 
Dunne et al. (2003). 
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FIG. 5. — The hatched region shows the evolution of the stellar mass density, as in Figure 3. The region hatched with dashed lines is p gas (z) = p*(0) + 
Pgas(0) — p*(z) (see text for details). The open points indicating total baryonic density at z = (Fukugita & Peebles 2004), z = 2.5 (Dunne et al. 2003), and 
z = 3 (Fukugita et al. 1998). The filled circles are the neutral gas mass density (H+He) in DLA systems. 



